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In a recent preprint James Owen Weatherall has attempted a simple local-deterministic model for 
the EPR-Bohm correlation and speculated about why his model fails when my counterexample 
to Bell's theorem succeeds. Here I bring out the physical, mathematical, and conceptual reasons 
- - - why his model fails. In particular, I demonstrate why no model based on a tensor representation 

m : of the rotation group SU(2) can reproduce the EPR-Bohm correlation. I demonstrate this by 

calculating the correlation explicitly between measurement results iS' = ± 1 and S§ = ±1 in a. local 
C " ^ ■ and deterministic model respecting the spinor representation of SU(2). I conclude by showing how 

' Weatherall's reading of my model is misguided, and bring out a number of misconceptions and 

, , ' unwarranted assumptions in his imitation of my model as it relates to the Bell-CHSH inequalities. 
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I. INTRODUCTION 



In a recent preprint [1| James Owen Weatherall has attempted a local and deterministic model for the EPR-Bohm 
' experiment, which he affirms to have been inspired by my work on Bell's theorem 

^ [3^ H 1 • Although his analysis 

^ _ mainly concerns his attempted reconstruction of my model and why the attempt fails, his preprint has been worded in 
j^JT)' a manner that has allowed some readers to embrace his discussion as a criticism of my work on Bell's theorem Q ■ 
^ ■ Here I demonstrate that (as he himself stresses to some extent) the analysis Weatherall presents in his preprint has 
Q nothing to do with my model, or with the physics and mathematics of the EPR-Bohm correlation QQ- In fact his 
' ^ ' analysis exhibits no understanding of how my local-realistic framework works, nor of the reasons why it explains the 
J>^^ origins of all quantum correlations I show this by first calculating the EPR-Bohm correlation in a successful 

■ local-deterministic model based on the spinor representation of SU(2), and then revealing a number of misconceptions 
Q-i' and unwarranted assumptions in Weatherall's reconstruction of my model as it relates to the Bell-CHSH inequalities. 
I conclude that, contrary to first impressions, Weatherall's thinly veiled criticism of my work is entirely vacuous. 



II. AN EXACT LOCAL-DETERMINISTIC MODEL FOR THE EPR-BOHM CORRELATION 

In order to bring out the erroneous assumptions in Weatherall's analysis, it would be convenient to assess it in the 
light of a successful model for the EPR-Bohm correlation. This will allow us to unveil his assumptions more easily. 



A. A Complete Specification of the Singlet State 

To this end, let us recall Bell's definition of a complete theory [9]. He considered a physical theory to be complete 
just in case its predictions for the EPR-Bohm experiment are dictated by local-deterministic functions of the form 

^(r. A) : IR3 X A — ^ 5° = {-1, +1}, (1) 

where IR"^ is the space of 3- vectors r, A is a space of the complete states A, and S'^ = { — 1, +1} is a unit 0-sphere. He 
then proved a mathematical theorem concluding that no pair of functions of this form can reproduce the correlation 
as strong as that predicted by quantum mechanics for the rotationally invariant singlet state: 



f (a, b) = lim 

n 3> 1 



n 

- V ^(a, A'^) i^(b, A* 

n ^ — ^ 



n 
fe = i 



^-a-b. (2) 
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As it stands, this conclusion of Bell is entirely correct and beyond dispute, provided we accept prescription ((T]) as 
codifying a complete specification of the singlet state. In his preprint, following Bell and CHSH Weatherall 
accepts dl]) as codifying a complete specification of the singlet state, whereas my work begins by recognizing that ([T]) 
does not, and cannot, codify a complete specification of the singlet state I have argued that Bell's prescription is 

based on an incorrect underpinning of both the EPR argument ^llj and the actual topological configurations involved 
in the EPR-Bohm experiments Q , even if we leave the physics and mathematics underlying the correlation aside. My 
argument is rather subtle and requires a clear understanding of what is meant by both a function in mathematics and 
the geometry and topology of a parallelized 3-sphere. But the bottom line of the argument is that, for any two-level 
system, the EPR criterion of completeness demands the correct measurement functions to be necessarily of the form 



± 1 = £/(r. A) : R^x A — ^ 5^ - SU(2), 



(3) 



with the simply- connected codomain of A) replacing the totally disconnected codomain S'^ assumed by Bell. 
Thus J2/(r, A) = ± 1 now represents a point of a parallelized 3-sphere, S^. As a function, it takes values from the 
domain JR^ x A and ends up belonging to the codomain S^. Consequently, any correlation between a pair of such 
results is a correlation between points of a parallelized 3-sphere. Unless based on a prescription of this precise form, 
any Bell-type analysis simply does not get off the ground, because without completeness there can be no theorem f^]- 

Here — which can be thought of as the configuration space of all possible rotations of a rotating body (including 
spinorial sign changes) — is defined as the set of all unit quaternions isomorphic to a unit parallelized 3-sphere: 



q{^, r) 



exp 



/3(r) 



q(^, r)||^ 



(4) 



where ^(r) is a bivector rotating about r € IR^ with the rotation angle ip in the range < tp < Air. Throughout this 
paper I will follow the concepts, notations, and terminology of geometric algebra Accordingly, ^(r) G 5^ C S*^ 

can be parameterized by a unit vector r = ri ei + r2 62 + 63 e ]R^ as 



/3(r) (/.r) 

= ri ( / • ei ) + r2 ( / • 62 ) + ra ( / • 63 ) 
= ri 62 A 63 + r2 ea A ei + ei A 62 



(5) 



with /3^(r) = —1. Here the trivector / := ei A 62 A 63 (which also squares to —1) represents a volume form of the 



physical space 12| [13| . Each configuration of the rotating body can thus be represented by a quaternion of the form 

q(-0, r) = cos ^ -f /3(r) sin ^ , (6) 
which in turn can always be decomposed as a product of two bivectors belonging to an <Z ^ SU(2), 

/3(r")/3(r' 



t 
2 



cos — -t- /3(r) sin 



± 
2 ' 



(7) 



with ip being its rotation angle from q(0, r) = 1. Note also that q('0, r) reduces to ± 1 as i/' — ^ 2kh for k = 0, 1, or 2. 

It is of paramount importance to note here that our topologically corrected prescription ^ does not alter the actual 
measurement results. For a given vector r and an initial state A, both operationally and mathematically we still have 



^(r. A) = + 1 or 



(8) 



as the image points of the function ^ (r. A) as demanded by Bell, but now the topology of its codomain has changed 
from a 0-sphere to a 3-sphere, with the latter embedded in IR^ in such a manner that the prescriptions ([1]) and ^ 
are operationally identical ^ [w] . On the other hand, without this topological correction it is impossible to provide a 
complete account of all possible measurement results in the sense specified by EPR. Thus the selection of the codomain 
^ IR* in prescription ([3]) is not a matter of choice but necessity. What is responsible for the EPR correlation is 
the topology of the set of all possible measurement results. But once the codomain of ^(r, A) is so corrected, the proof 
of Bell's theorem (as given in Refs. [9]) simply falls apart. Moreover, it turns out that the strength of the correlation 
for any physical system is entirely determined by the torsion within the codomain of the local functions £^{v. A). 

Returning to Weatherall's analysis, it should be clear now that, because it is based on a prescription other than ([3]), 
it too is a non-starter Q . More importantly, once we recognize that the only way of providing a complete account of all 
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possible measurement results for the singlet state is by means of prescription ([3]), the statistical procedure for analyzing 
the correlation must be consistently customized for the set of all unit quaternions, which Weatherall fails to do. 
Since this procedure can be appreciated more readily by studying the explicit construction of my model from my book 
and elsewhere, I now proceed to reproduce the model in some detail in the following subsections. Doing so will also 
dispel a persistent but gravely disingenuous charge that my model encounters "certain technical complications" (^j ^] . 



B. Construction of the Measurement Functions 

Once the measurement results are represented by functions of the form Q , it is easy to reproduce the EPR-Bohm 
correlation in a manifestly local, realistic, and deterministic manner. This is because a parallelized 3-sphere has quite 
a unique and distinctive topological structure [6] [7] . It is one of the only two parallelizable spheres with non- vanishing 
torsion — the other one, with variable torsion, being the 7-sphere. Once parallelized by a constant torsion, the 3-sphere 
remains closed under multiplication, and forms one of the only four possible normed division algebras Q Q ■ These 
are profound concepts underlying the very existence and strength of quantum correlations [3] . By ignoring them and 
dismissing them as irrelevant, Weatherall is ignoring the physics and mathematics of the quantum correlations. More 
importantly, because of the unique and distinctive topological characteristics of the 3-sphere Q, the measurement 
functions such as =2/(r. A) for the EPR-Bohm correlation have to be constructed in a very specific manner for any 
model to be successful To this end, let us begin with the following definition of the orientation of a vector space: 

Definition II. 1 An orientation of a finite dimensional vector space Vd is an equivalence class of ordered basis, 
say {/i, . . . , fd}, which determines the same orientation of Vd as the basis {/{, . . . , f'f} if f- — ujijfj holds with 
det{u)ij) > 0, and the opposite orientation ofVd as the basis {/{, . . . , /^} if f^ = t^ijfj holds with det(ujij) < 0. 

(Here repeated indices are summed over.) Thus each positive dimensional real vector space has precisely two possible 
orientations, which (rather suggestively) can be denoted as A = -fl or A = — 1. More generally an oriented smooth 
manifold such as consists of that manifold together with a choice of orientation for each of its tangent spaces. 

It is important to note that orientation of a manifold is a relative concept jl3 |. In particular, the orientation of a 
tangent space Vd of a manifold defined by the equivalence class of ordered basis such as {/i, . . . , fd} is meaningful 
only with respect to that defined by the equivalence class of ordered basis {/{, . . . , /^}, and vice versa. To be sure, 
we can certainly orient a manifold absolutely by choosing a set of ordered bases for all of its tangent spaces, but the 
resulting manifold can be said to be left or right oriented only with respect of another such set of ordered basis [14|. 

Now the natural configuration space for an EPR-Bohm type experiment is a unit parallelized 3-sphere, which can 
be embedded in IR** with a choice of orientation, say A = +1 or —1. This choice of orientation can be identified with 
the initial state of the particle pair in the singlet state with respect to the orientation of the detector basis as follows. 
We first characterize the embedding space IR'* by the graded basis 

{1, ii(A), L2(A), i3(A)}, (9) 

with A = ± 1 representing the two possible orientations of and the basis elements Lp{X) satisfying the algebra 

Lp{\)Ly{\) = -g^y - e^^pLpiX), (10) 

with an arbitrary metric g^^ on S^. Here the bivectors { Lp{\) } will represent the spin angular momenta of the 
particles, with /i = 1, 2, 3 and the repeated indices summed over. These momenta can be assumed to be detected by 
the detector bivectors, say { Dp }, with the corresponding detector basis {1, _D2, -D3 } satisfying the algebra 

DpDi, = -gpy - CpypDp (11) 

and related to the spin basis { 1, £i(A), L2{\), L-i{\) } as 



( ^ \ 












/ 1 






A 







Di 









A 




D2 


Vi3(A)/ 
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(12) 



Evidently, the determinant of this matrix works out to be det{u!ij) — A. Since A^ = -1-1 and is a 4 x 4 identity 
matrix, this relation can be more succinctly written as 



L^(A) = Ai^^ and i^^ = AL,,(A), 



(13) 
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or equivalently as 

{1, Li(A), L2(A), L^iX)} = {1, XDi, XD2, XD3} (14) 

and 

{1, D2, D3} = {1, Aii(A), AL2(A), XLsiX)}. (15) 

These relations reiterate the fact that orientation of any manifold is a relative concept. In particular, orientation of 5*^ 
defined by the spin basis { 1, L^{X) } is meaningful only with respect to that defined by the detector basis { 1, } 
with the orientation A = +1, and vice versa. Thus the spin basis are said to define the same orientation of as 
the detector basis if Lp(A = +1) = +0^^, and the spin basis are said to define the opposite orientation of S'^ as the 
detector basis if i^(A = —1) = — D^. Note also that the numbers 1 and L(r, A) are treated here on equal footing. 

We are now in a position to define the functions £/ (a, A) and 3i{h, A) as results of measurement interactions (or 
Clifford products) between detector bivectors — D(a) and +D(b) and spin bivectors L(a, A) and L(b, A) as follows: 

SU(2)^53 9±1 = ^(a, A'^) = -D(a)L(a, A'=) = {-a^D^}{a, L,(A'=)} = <'^; !^ (16) 



if A*^ = - 1 



and 



SU(2)^53 3±1 ^ ^(b, A'=) = +D(b)L(b, A'^) = {+6^D^}{5,L.(A'=)} = <j , j x" = -1 ^^^^ 



where the relative orientation A is now assumed to be a random variable, with 50/50 chance of being +1 or — 1 at 
the moment of creation of the singlet pair of spinning particles. In what follows, I will assume that the orientation of 
defined by the detector basis { 1, Di^ } has been fixed before hand 30- Thus the spin bivector { L^{X) } is a 
random bivector with its handedness determined relative to the detector bivector { a,y Di, }, by the relation 

L(a, A) = { i^(A) } = A { a, Z?, } = A D(a), (18) 

where, as a direct consequence of the algebra (fTO|) with (7^,^ = S^iy, the bivectors L(a, A) satisfy the following identity: 

L(a, A) L(a', A) = - a • a' - L(a x a'. A). (19) 

Using these relations the spin detection events (fT6|) and (fTT]) follow at once from the algebras defined in (fTOj) and (fTT|) . 

Evidently, the measurement results £/ (a, A) and ^(b. A) as defined above, in addition to being manifestly realistic, 
are strictly local and deterministically determined numbers. In fact, they are not even contextual. Alice's measurement 
result ^(a. A) — although it refers to a freely chosen direction a — depends only on the initial state A; and likewise, 
Bob's measurement result ^(b. A) — although it refers to a freely chosen direction b — depends only on the initial 
state A. Let us also not overlook the fact that, as binary numbers, ^(a. A) = ± 1 and ^(b. A) = ± 1 are still points 
of a parallelized 3-sphere. To confirm this, recall that a parallelized 3-sphere is a set of unit quaternions of the form 



q'(V', r. A) := |a^- cos | + L(r, A^) sin | | 



(20) 



and a measurement result such as ^(a. A) = ±1 is a limiting case of such a quaternion constituting the 3-sphere: 

9 ± 1 = A) = lim £/(a, a'. A) 

= lim {-D(a)L(a', A)} 

a' —> a 

= lim {(-/•a)(A/-a')} 

a' —> a 

= lim { A a • a' + A / • (a X a') } 

a'— > a 

= lim I A cos — + L(c, A) sin — 

V^2K7r \ 2 ^ ' ^ 2 

= hm {q(V', c. A)}. (21) 

-0— > 2K7r 

Here / = ei A 62 A 63 is the volume form, limit a' a is equivalent to the limit ip — 2ktt for k = 0, 1, or 2, 
tp = 2 77aa' is the rotation angle about the axis c := a x a'/|a x a'|, and r^aa' is the angle between a and a' Q. 
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C. A Crucial Lesson from Basic Statistics 



It is important to note that the variables ^(a, A) and ^(b, A) defined in equations (fT6| and (fT7| are generated 
with different bivectorial scales of dispersion (or different standard deviations) for each measurement direction a and 
b. Consequently, in statistical terms these variables are raw scores, as opposed to standard scores [T5l|. Recall that a 
standard score, z, indicates how many standard deviations an observation or datum is above or below the mean. If x 
is a raw (or unnormalized) score and x is its mean value, then the standard (or normalized) score, z(x), is defined by 



z(x) 



a(x) 



(22) 



where cr(x) is the standard deviation of x. A standard score thus represents the distance between a raw score and 
population mean in the units of standard deviation, and allows us to make comparisons of raw scores that come from 
very different sources In other words, the mean value of the standard score itself is always zero, with standard 

deviation unity. In terms of these concepts the correlation between raw scores x and y is defined as 



f(x, y) 



lim 



1 " 

- V (x'^- - x) (y'^ - y) 

77 L ^ 



n 

L fc = 1 



cr(x) cr(y) 



lim 

n> 1 



1 " 

- z(x'=) z(y'=) 



n 

L fe = 1 



(23) 
(24) 



It is vital to appreciate that covariance by itself — i.e., the numerator of equation (j23p by itself — does not provide the 
correct measure of association between the raw scores, not the least because it depends on different units and scales 
(or different scales of dispersion) that may have been used in the measurements of such scores. Therefore, to arrive 
at the correct measure of association between the raw scores one must either use equation (|23p . with the product of 
standard deviations in the denominator, or use covariance of the standardized variables, as in equation (|24p. 



Now, as discussed above, the random variables ^(a, A) and S§{h, A) are products of two factors — one random 
and another non-random. Within the variable i2^'(a. A) the bivector L(a, A) is a random factor — a function of the 
orientation A, whereas the bivector — D(a) is a non-random factor, independent of the orientation A : 



j^(a. A) 
and ^(b. A) 



-D(a)L(a, A) 
-D(b)L(b, A) 



(25) 
(26) 



Thus, as random variables, £^ (a. A) and ^(b. A) are generated with different standard deviations — i.e., different sizes 
of the typical error. More specifically, ^(a. A) is generated with the standard deviation — D(a), whereas S§ih, A) is 
generated with the standard deviation -f D(b). These two deviations can be calculated as follows. Since errors in the 
linear relations propagate linearly, the standard deviation a{s^) of ^(a. A) is equal to — D(a) times the standard 
deviation of L(a, A) [which I will denote as (j{A) — cr(La)], whereas the standard deviation (t(^) of !3§{h. A) is equal 
to -|-D(b) times the standard deviation of L(b, A) [which I will denote as cr(i?) = (T(Lb)]: 



and 



a{.s/) = -D(a)(T(A) 
a{.^) = +D(b)a(B). 



(27) 
(28) 



But since the bivector L(a, A) is normalized to unity, and since its mean value TO(La) vanishes on the account of A 
being a fair coin, its standard deviation is easy to calculate, and it turns out to be equal to unity: 



o(A) 



\ 



1 II 

fc=i 



A{sL, A'=) 



\ 



1 " 

-E l|L(a, A^) 



(29) 



where the last equality follows from the normalization of L(a, A). Similarly, it is easy to see that a{B) is also equal 
to 1. Consequently, the standard deviation of j2/(a. A) = ± 1 works out to be — D(a), and the standard deviation of 
.^(b. A) = ± 1 works out to be -f D(b). Putting these two results together, we arrive at the following standard scores 
corresponding to the raw scores and 



(30) 
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and 

^(b, A) - ^(bTX) +D(b)L(b, A) - 

= ^ Tn(b) = (^^^ 

where I have used identities such as — D(a)D(a) = +1. Needless to say, these standard scores are pure bivectors: 

SU(2) ^ S^D 3 L(a, A) = ± 1 about a e ]R^ (32) 
and SU(2) D 3 L(b, A) = ± 1 about b e (33) 



D. How Errors Propagate within a Parallelized 3-sphere 



As noted towards the end of subsection II A, one of several oversights in Weatherall's reading of my model concerns 
his failure to recognize the necessity of applying the correct statistical procedure for analyzing the correlation between 
the measurement results defined by (IT6l) and (ITTI) . This is surprising, because I pointed out this oversight to him in 
a private correspondence more than eighteen months ago. As we noted above, a parallelized 3-sphere is a set of unit 
quaternions. Each point of a parallelized 3-sphere is thus represented by a unit quaternion. As a result, the correct 
statistical procedure within my model must take into account how errors propagate in a 3-sphere of unit quaternions. 

Accordingly, let a probability density function P(q) : [0, 1] of random quaternions over be defined as: 

P(q) = , ' J ||q-m(q)||n 

/2.||a(q)||^ I 2||a(q)|r j 

where the square root of q = pp, p e S^, is defined as 

yq = ypp := ±pt(pp) = ±(ptp)p = ±p. (35) 

It is a matter of indifference whether the distribution of q e S*"^ so chosen happens to be normal or not. Here q is 
an arbitrary quaternion within 5*^ (A) of the form (1201) . which is a sum of a scalar and a bivector (treated on equal 
footing), with < ip < At: being the double-covering rotation angle about r-axis. The mean value of q is defined as 

1 " 

m(q) = - ^ q\ (36) 

fc = 1 

and the standard deviation of q is defined as 
CT[q(V', r, A)] 



i J2 {(iH^) - "i(q)} {q'=(27r-V') - m{q)}^ . (37) 
\ fe = 1 



Note that in this definition q(V') is coordinated by to rotate from to 2tt, whereas the conjugate q' {2tt — -0) is 
coordinated by ip to rotate from 27r to 0. Thus, for a given value of A, both and q^^(27r — ip) represent the same 
sense of rotation about r (either both represent clockwise rotations or both represent counterclockwise rotations). 
This is crucial for the calculation of standard deviation, for it is supposed to give the average rotational distance 
within from its mean, with the average being taken, not over rotational distances within a fixed orientation of S^, 
but over the changes in the orientation A of itself. Note also that, according to the definition ((20)) . q{ip) and its 
conjugate q^(V') satisfy the following relation: 

qt(27r-V) - -q(V^). (38) 

Consequently, the standard deviation of both qt(27r — tp) and — q(V') must, with certainty, give the same number: 

a[qt(2^-0)] = a[-q(V.)]. (39) 

It is easy to verify that definition ([37|) for the standard deviation of does indeed satisfy this requirement, at least 
when TO(q) = 0. What is more, from equation (j38|) we note that the quantity being averaged in the definition (|37p 
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is essentially — q q. This quantity is insensitive to spinorial sign changes such as q — — q, but transforms into the 
quantity — q^q^ under orientation changes such as A — >■ —A. By contrast, the quantity — qq^ would be insensitive 
to both spinorial sign changes as well as orientation changes. Thus tT[q(A)], as defined in (|37| . is designed to remain 
sensitive to orientation changes for correctly computing its averaging function on q(A) in the present context. 

Now, in order to evaluate ct(j2/) and cr(La), we can rewrite the quaternion (I20p rotating about r = a as a product 

qi^jj, a, A) = p{^P, a) L(a, A) (40) 

of a non-random, non-pure quaternion 

p(V,a) ■.= sm(t-) _ D(a) cos = exp |- D(a) f^)! (41) 



^2 J ' ' V 2 
and a random, unit bivector L(a, A) satisfying 



1 " 

- ^ L(a, A'=)L'f(a, A*^) = 1 . (42) 



n 

k = l 



Note that p{ip, a) reduces to the unit bivector ± D(a) for rotation angles ip = 0, i/j — 2tt, and ijj = An. Moreover, using 
the relations L(a, A) — AD(a) and D^(a) = — 1 it can be easily checked that the product in (^0]) is indeed equivalent 
to the quaternion defined in ((20|) for r = a. It is also easy to check that the non-random quaternion p{ip, a) satisfies 
the following relation with its conjugate: 

pt(27r- V, a) = p(V', a). (43) 

Consequently we have 

qt(27r - 7/., a. A) {p(27r - t/-, a) L(a, A)}^ = Lt(a, A) pt(27r - t/j, a) = V {a, A) p(V^, a) . (44) 
Thus, substituting for q(?/', a. A) and q1'(27r - ip, a. A) from Eqs. gOl) and (j44]) into Eq. ([371), together with 



"^(^a) = ^ E qa = P(^, a) ji L(a, A'^)) = p(V., a) jl a4 D(a) = 0, 

It \ I L I III/ I 

k = l (, k = l ) \ k = l ) 



(45) 



we have 

cr[q(V', a. A)] = 



n 

{p(^, a)L(a, A''0}{Lt(a, A'=)p(^,a)} 
\ fc = 1 



P(*, «) I ^ I'fa. ^')Lt(a, A') | p(V., a) 





= V P(i/', a)p(V', a) 

-±P(V', a). (46) 

Here I have used the normalization of L(a, A) as in (UU, and the last equality follows from the definition ([55)) . It can 
also be deduced from the polar form of the product 

p(?/', a) p(?/;, a) = cos (tt — i/i) — D(a) sin [tt — ij;) — exp { — D(a) (tt — ?/;)} . (47) 

The result for the standard deviation we have arrived at, namely 

a[q(V', a. A)] = ±p(V', a), (48) 

is valid for all possible rotation angles t/i between the detector bivector — D(a) and the spin bivector L(a, A). For the 
special cases when -0 = 0, tt, 27r, Stt, and 47r, it reduces to the following set of standard deviations: 

ct[ q(V' = 0, a. A)] = o-(£/) = ± D(a) 

Cr[q(-0 = TT, a. A)] = cr(La) = ± 1 

cr[ q^ijj = 27r, a, A)] = cr(i/) = ± D(a) 

= Stt, a, A)] — cr(La) = ± 1 

and f7[q(V' =47r, a. A)] = cr(i/) = ±D(a). (49) 
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To understand the physical significance of these results, let us first consider the special case when ip = tt. Then 

q(7/^ = ^, a, A) = +L(a, A), (50) 

which can be seen as such from the definition (|20l) above. Similarly, for the conjugate of q{ip = vr, a. A) we have 

qt(i^ = ^, a. A) = -L(a, A) = + V{a, A). (51) 

Moreover, we have m(La) = 0, since L(a, A) = + AD(a) with A = ± 1 being a fair coin. Substituting these results 
into definition (|37|) — together with tp = tt — we arrive at 



CT(La) = 



1 " 

L(a, Afe)Lt(a, A'^) = ±1, (52) 
\ fe = i 



since L(a, A) Lt(a, A) = 1. Similarly, we can consider the case when tp = Sn and again arrive at cr(La) = ± 1. 

Next, we consider the three remaining special cases, namely "0 = 0, 27r, or Air. These cases correspond to the 
measurement results, as defined, for example, in equation (|16p. To confirm this, recall from equation ()2ip that a 
measurement result such as j2/(a. A) = ±1 is a limiting case of the quaternion If we now rotate the bivector 

L(c, A) to L(a, A) as 

D(r)L(c, A)Dt(r) = L(a, A) (53) 
using some D(r), and multiply Eq. (I^TI) from the left by D(r) and from the right by D^(r), then we arrive at 

i^(a, A) = lim {D(r)q(7/;, c, A)Dt(r)} 

a'— V a 

= lim I A cos — + L(a, A) sin — 

>2K7r 2 2 

= lim { q(V', a. A) } 

= lim {p(0, a)L(a, A)}, (54) 

where p(V', a) is defined in equation (HIT) . The limit a' -> a is thus physically equivalent to the limit ip — >■ 2ktt for 
K = 0, 1, or 2. We therefore have the following relation between A) and a. A): 

q(V' = 2ktt, a. A) = ± D(a) L(a, A) = ± £/{&, A) , (55) 

and similarly between A) and q^(?A, a. A): 

qt(?A = 2ktt, a. A) = { ± D(a) L(a, A) }^ = ± £/\a, A) . (56) 

For example, for -0 = the definition (pil)) leads to 

q(0 = 0, a. A) = - D(a) L(a, A) = + ^/(a. A) . (57) 

This tells us that in the ip limit the quaternion q(V', a, A) reduces to the scalar point — D(a)L(a, A) of S^. 
Moreover, we have m{s^) — 0, since m(La) = as we saw above. On the other hand, from the definition (pUj) of 
q('0, a. A) we also have the following relation between the conjugate variables ja?'^(a. A) and q^(V' = 27r, a. A): 

qt(Vj = 27r, a. A) = + {D(a)L(a, A) }^ = +Lt(a, A)Dt(a) = -Lt(a, A)D(a) = -i3ft(a. A). (58) 

This tells us that in the — > 27r limit the quaternion q^{'^p^ a. A) reduces to the scalar point — Lt(a, A)D(a) of 
. Thus the case ip — does indeed correspond to the measurement events. The physical significance of the two 
remaining cases, namely ip — 211 and 47r, can be verified similarly, confirming the set of results listed in (|49p : 



cr(j^) = ±D(a) . (59) 
Substituting this and a{S§) = ±D(b) into Eqs. ((30)) and ((3T|) then immediately leads to the standard scores: 



±^(a. A) --^(a. A) ±D(a)L(a, A) - f ±D(a) , 



and PAh - ^) ~ ^) - ±D(b)L(b, A) - _ f ±D(b) ] _ 

This confirms the standard scores derived in the equations pOI) and (1311) of the previous subsection. 
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FIG. 1 Propagation of the 68% probability interval from a random bivector S to a scalar si within the parallelized 3-sphere. 



So far I have assumed that randomness in the measurement results i2^(a, A) and ^(b, A) originates entirely from the 
initial state A representing the orientation of the 3-sphere. In other words, I have assumed that the local interaction of 
the detector D(a) with the random spin L(a, A) does not introduce additional randomness in the measurement result 
^ (a, A). Any realistic interaction between D(a) and L(a, A), however, would inevitably introduce such a randomness, 
of purely local, experimental origin. We can model this randomness by introducing an additional random variable, 
say '"a € [0, 1], not dependent on A. Physically we can think of ra as an alignment parameter between the detector 
bivector D(a) and the spin bivector L(a, A), with ra = 1 representing the perfect alignment. Clearly, introduction of 
this additional random parameter will make all the bivectors and quaternions unnormalized, and the corresponding 
probability density function ([M)) would then represent a Gaussian distribution — provided we also assume that the 
orientation A = ± 1 of S*"^ itself is distributed non-uniformly between its values +1 and —1. Moreover, although the 
measurement results would then fall within the range — 1 < ^(a. A) < +1, their mean value would be zero for a 
uniformly distributed A, since the mean value of the product of the independent random variables ra and A would 
then be the product of their mean values: 

tnira.X) = m(ra)m(A). (62) 

More importantly, the standard scores computed above in equations ([50)1 and (pT|) would not be affected by this more 
realistic random process ra A — at least for the special case of uniformly distributed A — because they involve the ratios 
of the corresponding raw scores and standard deviations centered about the mean values mirg. A) = = m{r\y A): 



±^(a. A) - ^(a. A) ±raD(a)L(a, A) - \ ±raD(a) ) ^ 

^) = ^) = ^) = 1 .K) I ^) = ^) 

, ±/^(b. A) - ^(b. A) ± rbD(b)L(b, A) - r±rbD(b)1 „ .^ 

Let us now try to understand the propagation of error within from this physically more realistic perspective. 
To this end, let the random variable q('!/') a, A) € S^{X) be such that the measurement results ^(a. A) S [—1, +1] 
remain as before, but the bivectors L(a, A) are subject to a random process raA such that S(a, A, ra) = raL(a, A) 
with ra G [0, 1]. Then the mean value m(S) and standard deviation (t(S) of S would be a bivector and a scalar: 

m(S) = a bivector 
and (t(S) = a scalar. (65) 

If we now take the detector bivector to be D(a) = / • a as before, then the measurement results can be identified 
as — 1 < ^ = DS < +1 so that m{£/) > 0. Since D is a non-random bivector, errors generated within jz/ by the 
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random process A would stem entirely from the random bivector S, and propagate linearly. In other words, the 
standard deviations within the random number £/ due to the random process rg A would be given by 

ct(^) = D(7(S). (66) 

But since cr(S), as we noted, is a scalar, the typical error a{£/) generated within £/ due to the random process Ta A 
is a bivector. The standardized variable (which must be used to compare the raw scores with other raw scores ^ ) 
is thus also a bivector: A :— /cr{£/) — scalar x S. 

As straightforward as it is, the above conclusion may seem unusual. It is important to recall, however, that in 
geometric algebra both scalars and bivectors are treated on equal footing [l2|[i3- They both behave as real- valued 
c-numbers, albeit of different grades. To appreciate the consistency and naturalness of the above conclusion, let 

= /(S) = DS (67) 

be a continuous function generated by the geometric product of two bivectors D(a) and S(a, A, Ta) as before. The 
natural question then is: How does a typical error in S governed by the probability density (j34p — which can be 
represented by the 68% probability interval 

[m(S) - o-(S), m(S) + cr(S)] (68) 

as shown in the Fig.[T] — propagate from the random bivector S to the random scalar through the function 
/(S) = D S? To answer this question we note that the two end points of the interval ((68)) represent two points, 
say and q+, of the 3-sphere, which is a Riemannian manifold. The geometro- algebraic distance between the points 
and q+ can therefore be defined, say, as 

d{cr, q+) ^ (q" - q+) x sign(q- - q+) . (69) 

Moreover, from definition (j67p of and a first-order Taylor expansion of the function /(S) about the point S = to(S) 
we obtain 



^ = f{m{S)) + II 



(S-m(S)) + ... (70) 

S= m(S) 



Now it is evident that the slope df /dS = T) of this line is a constant. Therefore the mean m(j2/) and the standard 
deviation (t{£/) of the distribution of jz/ can be obtained by setting S = m(S) and S = cr(S): 

m(£/) = /(m(S)) = Dm(S) = a scalar (71) 
df 

and cr(^) = cr(S) = D a(S) = a bivector. (72) 
oS 

The probability distribution of jz/ is thus represented by the 68% interval 

[m{£/) - a{£/), to(£/) + ct(j2/) ] . (73) 

If we now set Ta = 1 and thereby assume that S is in fact the unit bivector L with a vanishing mean, then we have 
m{jz/) — and a{£/) = ± D, as in equation ([5^ above. 

Finally, it is instructive to note that, geometrically, the propagation of error within is equivalent to a simple 
change in the perspective (cf. Fig.[T]): 

bivcctor scalar scalar bivcctor 

3 to(S) ± cr(S) ^^^^ > m{£/) ± a{£/) . (74) 

V V 

quaternion quaternion 



In particular, the probability density of the scalar js/ over corresponding to interval ([73|l is equivalent to that of the 
bivector S over corresponding to interval (|68p . With this, we are now ready to derive the EPR-Bohm correlations. 
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E. Derivation of Pair Correlations Among the Points of 



We begin by noting that, according to my model, EPR-Bohm correlations are correlations among the points of a 
parallelized 3-sphere Q. Now, since we have assumed that initially there was 50/50 chance between the right-handed 
and left-handed orientations of the 3-sphere {i.e., equal chance between the initial states A = -I- 1 and A = — 1), the 
expectation values of the raw scores £/ (a. A) and ^(b. A) vanish identically. On the other hand, a.s discussed above, 
the correlation between these raw scores (or their first product moment coefficient d la Pearson [l5|) can be obtained 
only by computing the covariance between the corresponding standardized variables ^(a. A) and B(h, A), which gives 



f(a, b) 



lim 

n» 1 



lim 

n> 1 



n 



E 

fc= 1 



A(a, A'=)B(b, A*^) 



' Qfif bi, 



lim 

n> 1 



lim 

0, 



1 



n 



L fc = 1 

n 

n ^ 

k = l 



(75) 



where I have used algebra defined in pH)) and the relation (IT51) . Consequently, as explained in the paragraph just 
below Eq. ([24]), when the raw scores — ±1 and ^ = ± 1 are compared, their product moment will inevitably yield 



£ (a, b) lim 

n > 1 



n 

E 

k = l 



X'') ^(b, A*^) 



(76) 



since the correlation between the raw scores and is equal to covariance between the standard scores A and B. 

So far in this section we have put no restrictions on the metric tensor, which, in the normal coordinates centered 
at a point of would be of the form 



.x^x^ + 0{\x\^) 



(77) 



In other words, the algebra PU)) we have used in the derivation of correlation (1751) is a general Clifford algebra, with 
no restrictions placed on the quadratic form [l6j . On the other hand, if the codomain of the measurement functions 
^(a, A) is taken to be a parallelized 3-sphere, then the above metric tensor specializes to the Euclidean metric Sfj,^, 
because the Riemann curvature tensor of a parallelized 3-sphere vanishes, inducing a non- vanishing torsion This 
case corresponds to the geometry of the group SU(2) and specializes the correlation ([75)) to exhibit maximum strength: 



£{a, b) = -g^^a^b^ 



^ jJsV ^fJ, 



cos ?7ab ■ 



(78) 



which in turn manifests the sensitivity of £/ (a. A) and ^(b, A) to spinorial sign changes. To appreciate the significance 
of these changes @ , recall from subsection II A that a parallelized 3-sphere is a set of unit quaternions of the form 



q(i^, r) = cos ^ + /3(r) sin ^ , (79) 

with ip being the rotation angle. It is easy to check that q{ip, r) respects the following rotational symmetries: 

q(V' + 2K7r, r) = - q(i/), r) for k = 1, 3, 5, 7, . . . (80) 
q(V' -(- 4K7r, r) = -f q(7/;, r) for k = 0, 1, 2, 3, . . . (81) 

Thus q^{ip, r) correctly represents the state of a body that returns to itself only after even multiples of a 27r rotation. 

It is very important to appreciate that the strong correlation derived in ([75]) are correlation among the points of a 
parallelized 3-sphere, S^, taken as the codomain of the measurement functions s^{r, A). Thus the strength and the 
very existence of the EPR-Bohm correlation (or of any correlation for that matter) stem entirely from the topological 
properties of the codomain of the measurement functions si/ (r, A). Had we chosen the codomain of s/ (r, A) to be any 
manifold other than a parallelized 3-sphere, the resulting correlation would not have been as strong as — cosr/ab- 
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F. Derivation of Upper Bound Exceeding the Bell-CHSH Bound 



Returning to the expectation value (|76|) in its most general form we can now proceed to derive the Bell-CHSH-type 
bound on possible correlations (3|[i3]- To this end, consider four observation b, and b', for the standard 

EPR-Bohm experiment. Then the corresponding CHSH string of expectation values [3|, namely the coefficient 



5(a, a', b, b') := 5(a, b) + £(a, b') + £(a', b) - f (a', b') , 



(82) 



would be bounded by the constant 2-^2, as discovered by Tsirel'son within the setting of Clifford algebra applied to 
quantum mechanics in general 0][l3- Here each of the joint expectation values of the raw scores ^(a, A) = ± 1 and 
^{h, A) = ± 1 are defined as 



£(si, h) = lim 



n 

- £/{a, A'^) ^(b, A'^ 



n 

L fe = 1 



with the binary numbers such as ^ (a, A) defined by the limit 

3±1 = £/{a, A) = lim { q(V', a, A) } = - D(a) L(a, A). 



(83) 



(84) 



Thus £^ (a, A) and 3§{h, A) are points of a parallelized 3-sphere and £(a, b) evaluated in ((83|) gives correlation between 
such points of the 3-sphere [Sj. The correct value of the correlation, however, cannot be obtained without appreciating 
the fact that the number j2/(a. A) = ± 1 is defined as a product of a A-m-dependent constant, namely — D(a), and 
a A-dependent random variable, namely L(a, A). Thus the correct value of the correlation is obtained by calculating 
the covariance of the corresponding standardized variables 



A^{\) = A(a, A) = L(a, A) 
and Bb(A) = B{h, A) = L(b, A) , 



(85) 
(86) 



as we discussed just below equation ([24| . In other words, correlation between the raw scores (a, A) and ,^(b. A) is 
obtained by calculated the covariance between the standard scores A(a, A) and -B(b, A), as in equation ((75|) above: 



f (a, b) — lim 

n > 1 



1 " 

- ^ A(a, A^) B(b, A^) 



(87) 



The correlation between the raw scores is thus necessarily equal to the covariance between the standard scores: 



£(a, b) — lim 

n » 1 


- n 

- y ^(a, A^-) ^(b. A'^) 
. fc = i 


~ lim 

n» 1 


- n 

- y A(a, A'^) B(b, A'^) 
. fc = 1 





Using this identity we can now rewrite the CHSH string of expectation values (1821) in two equivalent expressions, 



iS(a, a', b, b') = lim 

n » 1 



n 

-y ^a(A^)^b(A^-) 



k = l 



lim 

n> 1 



1 " 

- y ^a(A'=)^b'(A'=) 



k = l 



lim 

n»l 



n 

L fe = i 



-J n 1 ^ 

-y ^a'(A'=)^b(A'^) - hm^ -y ^a'(A'=)^b'(A'' 



k = l 



(89) 



and 



iS(a, a', b, b') = lim 

n > 1 



- y Aa(A^-)i?b(A'^) + lim - y A,(A''li?b'(A'-) 



fe= 1 



k = \ 



lim 



n 

-y A^,{\^)B^{\^)) 



k = \ 



lim 

n» 1 



1 " 

-y A.,(A'=)i?b'(A^) 



(90) 



Our goal now is to find the upper bound on these strings of expectation values. To this end, we first note that the four 
pairs of measurement results occurring in the above expressions do not all occur at the same time. Let us, however. 
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conform to the usual assumption of counterfactual definiteness and pretend that they do occur at the same time, at 
least counterfactually, with equal distribution. This assumption allows us to simplify the above expressions as 



S{a, a', b, b') ~ lim 

n » 1 

and 



S{a, a', b, b') = lim 

n > 1 



n 

- V {<(A'=)^b(A'=) + =<(A'=)^b'(A^) + ^a'(A'=)^b(A'=) - <'(A'=)^b'(A'=)} 

n ^ — ^ ^ ^ 



k = l 



1 " 

- {^a(A^-)i?b(A'=) + Aa(A'=)i?b'(A'=) + ^a'(A")i?b(A'=) - ( A'^ ) Bb' ( A" ) } 



fc = l 



(91) 



(92) 



The obvious question now is: Which of these two expressions should we evaluate to obtain the correct bound on 
S{a, a', b, b')? Clearly, in view of the identity both expressions would give one and the same answer Thus it 
should not matter which of the two expressions we use to evaluate the bound. But it is also clear from the discussion 
in subsections II C and II D that the correct bound on the expression |9T|) involving the raw scores £/ and ^ can only 
be obtained by evaluating the expression (|92l) involving the standard scores A and B. Stated differently, if we tried to 
obtain the bound on S{a, a', b, b') by disregarding how the measurement results have been generated in the model 
statistically, then we would end up getting a wrong answer. By following the Bell-CHSH reasoning blindly Weatherall 
ends up making such a mistake. In the end S{a, a', b, b') is a functional of a random variable, and as such proper 
statistical procedure tailored to my model must be employed for its correct evaluation. This is an important point and 
the reader is urged to review the discussions in subsections II C and II D once again to appreciate its full significance. 

With these remarks in mind we proceed to obtain the upper bound on S{a, a', b, b') by evaluating the expression 
([5^ as follows. Since the standard scores v4a(A) — L(a, A) and -Bb(A) = L(b, A) appearing in this expression represent 
two independent equatorial points of the 3-sphere, we can take them to belong to two disconnected "sections" of 
(i.e., two disconnected 2-spheres within S^), satisfying 

[^r(A), Br'(A)] =0 V r and r' G 1R^ (93) 

which is equivalent to anticipating a null outcome along the direction r x r' exclusive to both r and r'. If we now 
square the integrand of equation (|92[) . use the above commutation relations, and use the fact that all bivectors square 
to —1, then the absolute value of S{a, a', b, b') leads to the following variance inequality [SJ: 

|5(a, a', b, b')| = \£{a, b) + £{a, b') + £{a', b) - £{a', b')| 



. lim 



E 

k = l 



{4 + 4raa'(Afc)rb'b(A'=)} 



where the classical commutators 



and 



raa'(A) :=-[Aa(A), Aa'(A)] = 



v4a 



'(A) 



(94) 



(95) 



rb'b(A) := -[Sb'(A), Bb(A)] = -Bb'xb(A) 



(96) 



are the geometric measures of the torsion within S^. Thus, it is the non- vanishing torsion T within the parallelized 
3-sphere — the parallelizing torsion which makes its Riemann curvature tensor vanish — that is ultimately responsible 
for the strong quantum correlation 3] 7]. We can see this at once from Eq. (j94| by setting T — 0, and in more detail 
as follows: Using definitions ([55]) and ([55)) for ^a(A) and -Bb(A) and making a repeated use of the bivector identity 

L(a, A) L(a', A) = - a • a' - L(a x a'. A) (97) 

specialized for the metric g^,^ = (5^^ on S^, the above inequality for 5(a, a', b, b') can be further simplified to 



|5(a, a', b, b')| s;; 




4 - 4 (a X a') • (b' X b) - 4 lim 

n > 1 



\ 



s; 2 a/ 1 - (a X a') • (b' X b) - 



n 

n ^ 

L fc = 1 



D(z) 



(98) 
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where z = (a x a') x (b' x b), and — as before — I have used the relation (|T8)) between L(z, A) and D(z) from subsection 
II B. Finally, by noticing that (a x a') • (b' x b) is bounded by trigonometry as 



-Is; (a X a') • (b' x b) ^ +1 , 
the above inequality can be reduced to the form 

- 2%/2 s$ f(a, b) + f(a, b') + £:(a', b) - £(a', b') s: +2\/2, 



(99) 
(100) 



which exhibits an extended upper bound on possible correlations. Thus, when in an EPR-Bohm experiment raw scores 
j2/ = ±1 and ^ — ±1 are compared by coincidence counts [l^ , the normalized expectation value of their product 



f (a, b) 



"C++ (a. 


b) 


h C_ 


-(a, b) - 


- c+ 


-(a, b) 


- C_+(a, b)" 


"C++ (a, 


b) - 


h C_ 


-(a, b) ^ 


h c+ 


-(a, b) 


+ C_+(a, b) 



(101) 



is predicted by my model to respect, not the Bell-CHSH upper bound 2, but the Tsirel'son upper bound 2\/2, where 
C-i (a, b) etc. represent the number of joint occurrences of detections + 1 along a and — 1 along b etc. 

This completes the presentation of my local, realistic, and deterministic model for the EPR-Bohm correlation. 



III. PHYSICAL, MATHEMATICAL, AND CONCEPTUAL FALLACIES OF WEATHERALL'S MODEL 



1. What is wrong with Weatherall's measurement ansatz? 



With the successful model firmly in place, we are now in a position to understand why Weatherall's model fails. 
To begin with, his model is based on a different representation of the rotation group. It is in fact based, not on 
the spinorial rotation group SU(2), but something akin to its tensorial cousin SO (3), which is a group of all proper 
rotations in IR"^, insensitive to spinorial sign changes. In fact, Weatherall takes a rather odd space, namely IR"^ A IR'^, 
for the codomain of his measurement functions ^(r. A), and then introduces another projection map to arrive at the 
measurement results {—1, +!}• Compared to my measurement functions (jl6l) obtained through a smooth limiting 
process (HJ), his two-step measurement process is rather artificial. Moreover, since his is not a simply-connected, 
parallelized codomain such as 5''^ that remains closed under multiplication, it cannot possibly satisfy the completeness 
criterion of EPR 3 • It is therefore not surprising why Weatherall is unable to find strong correlation among its 
points. Moreover, at the end of his two-step process his measurement results {—1, -f 1} are no longer the image points 
within the codomain IR"^ A IR'^ of his measurement functions. This is in sharp contrast with the situation in my 
model, where my measmenient results {—1, -f 1} remain very much a part of the codomain of my measurement 
functions. The reason why this comes about naturally within my model is because S^, the set of unit quaternions, is 
a simply-connected surface embedded in IR* that is equipped with a graded basis made of both scalars and bivectors: 

{1, 62 A 63, 63 A ei, ei A 63} . (102) 

Thus the scalars {—1, +1} are as much a part of as the bivectors L(r, A) are, regulated by these unified basis [3]. 
I am tempted to quip: What Nature has joined together, let no man put asunder. By contrast Weatherall's codomain 
is disconnected between the space IR'^ A IR'^ of "bivectors" and the set { — 1, -f 1} of scalars. His image points can thus 
be at best either bivectors or scalars, but not both. It is a disjoint world, more like the world of quantum mechanics. 

Let us, however, be more charitable to Weatherall. Let us grant him the codomain of his measurement functions to 
be the connected real projective space IRP'^, which is homeomorphic to the rotation group S0(3). After all, he does 
mention the Lie algebra so(3) in one of his footnotes. So let us grant him the smooth one-step measurement process 

± 1 = ^(r. A) : IR^ X A ^ IRP'^ - S0(3) (103) 

to reach the image subset {—1, +1}. This smooth map is well-defined within my model, since IRP"^ is simply the set 
of unit quaternions (cf. Eq. (|4])) with each point identified with its antipodal point [6,]. The measurement results 
± 1 e IRP'^ ^ SO (3) are thus limiting points of a quaternion, just as in equation (PT|) . The map that takes us from 
to IRP'^ can now be used to project the metric Sf^i, on onto IRP'^ to obtain the following induced metric on IRP"^: 

{- 1 + I 7?ab if < 77ab < TT 
(104) 
+ 3 - I ??ab if TT < 7?ab < 27r . 
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FIG. 2 Local-realistic correlations among the points of a parallelized 3-sphere can be stronger-than-classical but not quantum. 



Further details of how this metric is obtained from the metric on can be found in section III of Ref. [61] . The two 
metrics and J^ti/ thus provide relative measures of geodesic distances on the manifolds and IRP'^, respectively. 
Substituting the metric on HP^ into equation ([75)1 the correlation between the points of SO (3) then works out to be 

{- 1 + f ??ab if < ryab < TT 
(105) 
+ 3 - I ?yab if TT < ?7ab < 27r . 

The two sets of correlations, (|78|) and ()105p . are compared in Fig.[2j The general correlation function £{a, b) derived 
in equation (1751) can thus serve to distinguish the geodesic distances ^(a, b) on the groups SU(2) and S0(3) 



2. Why did we lose the strong correlation for S0(3)? 



It is crucial to appreciate that even when we choose S0(3) as the codomain of the function ^/{r, A) the correct 
statistical procedure that must be followed is the one described in subsections II C and II D above. This is because, 
as Weatherall himself notes, the Lie algebras of SU(2) and S0(3) are isomorphic to each other. In other words, the 
local algebraic or tangent space structures on SU(2) and S0(3) are identical, but not their metrical structures in the 
sense of geodesic distances. Thus the above statistical procedure, tailored to the graded basis (|102p . leading up to the 
general expression (|76p for correlations and beyond, is equally inevitable in the case of S0(3). Comparing the two 
sets of correlations resulting from this procedure — one for the prescription ([3]) and other for the prescription (jl03l) — it 
is then easy to see why we have lost the strong correlation in the second case. We started out with as a codomain 
of £/{r, A) and then, for the case of S0(3), we identified each point of with its antipodal point. But in doing so 
we lost the following spinorial rotation symmetry satisfied by q('0, r), as described in equations (|80p and (|8T|) above: 

q(V' + 2K7r, r) = - q(7/>, r) for k = 1, 3, 5, 7, . . . (106) 

In other words, by identifying the antipodal points of we lost the sensitivity to spinorial sign changes. As a result, 
q(i/', r) now represents the state of a rotating body that returns to itself after any and all multiples of 2tt rotation: 

q(i/' + 2ktt, r) = + q(V', r) for any k = 0, 1, 2, 3, . . . (107) 

This is the real reason why we lost the strong correlation for the SO (3) case. The reason Weatherall has argued for is 
an artifact of his bad choice of measurement functions. It stems from a failure to appreciate the unified nature of the 
graded basis (jl02l) and the associated fact that the scalars {—1, +1} and the bivectors L(r, A) occur as image points 
within the same codomain in my model. Thus the loss of correlation has nothing to do with the fact that ultimately 
the measurement functions must map to the image subset {—1, +1}. They manifestly do in my model [cf. Eqs. (jl6l) . 
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([TT)) . and (|2T]) ]. The raw scalar numbers £/ = ±1 and 3S = ±1 mapped to the image subset {—1, +1} — according to 
my model — are indeed the numbers used by Alice and Bob for calculating the correlation in the usual manner. And 
when, at the end of their experiment, they evaluate the statistical quantity 5(a, a', b, b') involving these numbers as 



iS(a, a', b, b') = lim 

n > 1 



n 

- V {<(A'=)^b(A'=) + ^4(A")^b'(A") + ^a'(A^')^b(A'^) - i/a'(A'=)^b'(A")} 



, (108) 



they will inevitably find that it exceeds the bound of 2 and extends to the bound of 2^/2. This conclusion may seem 
odd from the perspective of Bell-type reasoning, but the evidence presented for it in subsection II F is incontrovertible. 



IV. CONCLUDING REMARKS 



In addition to the main issue discussed above, it is instructive to reflect on the broader reasons why Weatherall's 
model fails. We can in fact identify at least six erroneous steps which engender the failure of his model from the start: 

1. Choice of incomplete codomain of the measurement functions ■-B/{r, A) (although with correct image points ±1). 

2. Neglect of putting scalars and bivectors on equal footing within a single, comprehensive, real number system. 

3. Failure to implement a spinor representation of SU(2) by recognizing the significance of spinorial sign changes. 

4. Lack of appreciation of the role played by the parallelizing torsion within S'^ for the existence and strength of 
strong correlations (or not recognizing the discipline of parallelization as the true source of strong correlations) . 

5. Failure to appreciate how errors propagate within S^, when taken as a codomain of the measurement functions. 

6. Neglect of the correct statistical procedure in the derivations of both the correlation and the upper bound 2\/2. 

Although interconnected, any one of these reasons is sufficient for the failure of Weatherall's model. Recognizing this, 
I must conclude that, contrary to first impressions, Weatherall's thinly veiled criticism of my work is entirely vacuous. 
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